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Abstract 



4^ ■ We present new imaging at 12.81 and 11.7 fim of the central ~40"x30" (~0.7x0.5 kpc) of the starburst galaxy 

M82. The observations were carried out with the COMICS mid-infrared (mid-IR) imager on the 8.2 m Subaru 
Q , telescope, and are diffraction-Hmited at an angular resolution of <0f'4. The images show extensive diffuse structures, 

■ including a 7''-long linear chimney-like feature and another resembling the edges of a ruptured bubble. This is the 
c/) I clearest view to date of the base of the kpc- scale dusty wind known in this galaxy. These structures do not extrapolate 

to a single central point, implying multiple ejection sites for the dust. In general, the distribution of dust probed in 
the mid-IR anticorrelates with the locations of massive star clusters that appear in the near-infrared. The 10-21 jim 
mid-IR emission, spatially-integrated over the field of view, may be represented by hot dust with temperature of 
^ ■ ~160 K. Most discrete sources are found to have extended morphologies. Several radio Hll regions are identified 

, for the first time in the mid-IR. The only potential radio supernova remnant to have a mid-IR counterpart is a source 

■ which has previously also been suggested to be a weak active galactic nucleus. This source has an X-ray counterpart 
in Chandra data which appears prominently above 3 keV and is best described as a hot (~2.6 keV) absorbed thermal 
plasma with a 6.7 keV Fe K emission line, in addition to a weaker and cooler thermal component. The mid-IR 
detection is consistent with the presence of strong [Nell] A 12.81/im line emission. The broad-band source properties 
are complex, but the X-ray spectra do not support the active galactic nucleus hypothesis. We discuss possible 
interpretations regarding the nature of this source. 

Key words: galaxies: starbursti — infrared: galaxies2 — techniques: high angular resolutions 
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I 1. Introduction 2008), and one new supernova (SN) is produced every ^3 years 

^ ■ (e.g. Rieke et al., 1980; Jones & Rodriguez-Espinosa, 1984). 

The galaxy M82 (NGC 3034) hosts the nearest and best ex- The energy output of the super star clusters hosting the SNe is 

ample of an ongoing massive starburst, making it an excellent thermalised and drives a large scale superwind along the galac- 

target for detailed studies at all wavelengths. The galaxy is tic minor axis (Heckman et al., 1990) which can be observed 

thought to have undergone an interaction event with its neigh- in detail because of the favorable edge-on inclination of the 

bor M81 about 10^ yr ago (Gottesman & Weliachew, 1977), source. 

triggering a massive nuclear starburst about 5x 10^ years ago In the infrared, M82 has been extensively studied with all 

(Rieke et al., 1980). There is also evidence for several other space missions. Its infrared luminosity is measured to be 

star formation episodes, both older and younger (de Grijs et al., 5 x 10^^ L©, and it shows prodigious dusty outflows and poly- 

2001; Forster Schreiber et al., 2003a). Around 40 supernova cyclic aromatic hydrocarbon (PAH) grains in the mid-infrared 

remnants (SNRs) have been identified in the core (Fenech et al., (mid-IR) extending on kpc-scales (e.g. Helou & Walker, 1988; 



Sturm et al., 2000; Forster Schreiber et al., 2003b; Engelbracht 

Based in part on data collected at Subaru Telescope, which is operated by ^ , n,r\r\^ j ^ ^ n,n,^r\ r> ^ ^ ^ n>f\^f\\ 

the National Astronomical Observatory of Japan. ^t al., 2006; Kaneda et al., 2010; Roussel et al., 2010). There 

are very few high- spatial-resolution mid-IR studies of the core 
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itself on scales of order 100 pc, because of size limitations of 
space missions. With their large primary mirrors, ground-based 
telescopes can achieve the best spatial resolution currently pos- 
sible. In the N (8-13 /im) band atmospheric window, the high- 
est resolution studies thus far are the works of Telesco & Gezari 
(1992, hereafter TG92) and Achtermann & Lacy (1995, here- 
after AL95), with nominal resolutions of li'l and 2'', respec- 
tively. 

In this work, we present the first sub-arcsec mid-IR N-band 
images of the core of M82. The galaxy was observed as an 
extension of our recent work of Seyfert galaxies (Gandhi et al., 
2009), as part of a study to understand the mid-IR emission 
of galaxies at high angular resolution. The observations were 
carried out at the 8.2 m Subaru telescope. At wavelengths of 
1 1 .7 and 12.8 1 /im, our imaging is diffraction-limited at < Of'4. 
These new images provide the sharpest mid-IR view of struc- 
tures at the base of the superwind, and allow an extensive multi- 
wavelength comparison of individual sources. Several Hll re- 
gions are newly identified in the mid-IR. We also discuss the 
nature of a putative active galactic nucleus candidate. Using a 
mid-IR detection and new Chandra data, we rule out the AGN 
hypothesis and discuss other possibilities, include supernova 
remnant ionization and emission from a starburst. 

Distance estimates to M82 have ranged over 3.2-5.2 Mpc 
(e.g. Burbidge et al., 1964; Tully, 1988; Sakai & Madore, 
1999). Some of the latest measurements suggest a distance 
at the lower end of this range based upon accurate determina- 
tion of the tip of the red giant branch magnitude (Karachentsev 
et al., 2004; Dalcanton et al., 2009). We adopt a value of 3.53 
Mpc herein, resulting in a physical scale of 17.1 pc per arcsec. 
Our imaging resolution limit corresponds to ^6.1 and 6.7 pc at 
11.7 /im and 12.8 /im, respectively. 

2. Observations 

Observations were carried out at Subaru on the night start- 
ing 2009 May 04, under generally clear weather conditions. 
The source was observed at the beginning of the night close to 
meridian, i.e. near its maximum altitude of 40°. The Cooled 
Mid-Infrared Camera and Spectrometer, or COMICS (Kataza 
et al., 2000) was used in imaging mode with standard chop- 
ping and nodding off- source. Due to the extended nature of the 
emission from M82, a relatively large chop throw of 30'' was 
used so that the sky position lay completely outside a single 
field-of-view of the detector. Technical difficulties prevented 
us from using larger throws. A North-South chop direction was 
adopted, because this is approximately perpendicular to the ap- 
parent major axis of the galaxy. Integration times of 60 s were 
used at each chop position. Subsequently, the telescope was 
nodded to a position at an offset of 1 arcmin to the North, and 
the above chopping series was replicated, this time entirely on 
sky. 

The above sequence was repeated 10 times, resulting in a to- 
tal on-source exposure of 600 s per filter and a total telescope 
time of about 50 min per filter after accounting for observing 
efficiency. Imaging was obtained in the Nell (narrow-band) 
and the Nl 1.7 filters, with central wavelengths (and full widths 
at half maximum) of 12.81 (0.2) /im and 11.7 (1.0) /im, re- 



spectively . The median airmass of the source was 1.56 and 
1.65, respectively. 

The Cohen standard star PI2 UMa (HD 73108; Cohen et al. 
1999) was observed for photometric calibration just before the 
target at an airmass of ^IA4. On-chip chopping with a throw 
of 10'' was used, and two nod observations were obtained with 
exposure time 10 s each. Seeing was measured in this obser- 
vation, as well as in a standard star observation of HD 108381 
(with identical observational setup to PI2 UMa; in the N11.7 
filter only) carried out approximately 1 h after the target ob- 
servation was completed. The maximum systematic error in 
the photometry due to the absolute calibration of the standard 
star is ^3%, and this error is added quadratically to all target 
photometric errors. 

3. Data reduction 

Data reduction was carried out using routines provided by 
the COMICS team (q_series package^), and following the 
procedures recommended in the COMICS Data Reduction 
Manual v.2.1.1. Some additional image manipulation was car- 
ried out in IRAF^. The COMICS data flow returns two sets 
of images for each observation: 1) a 'COMA' cube of data 
frames including all images obtained in one observing block; 
2) a single chop- subtracted 'COMQ' image, co-added over all 
the frames in that block. 

The recommended flat fielding procedure for imaging is a 
'self-sky-flat', in which the 'off beam chop positions are used 
for flat fielding the 'on' beam ones, and vice- versa. This is 
possible because the detector is illuminated quite uniformly by 
sky background, much brighter than the imaged sources. To 
create these flats, dark frames were first subtracted from each 
cube of data. For the standard star reduction, the on and the off 
chop beam data cubes were averaged to create a separate mean 
image for each of the two beams (hereafter, called 'mean beam 
images'). Low spatial frequency trends across the detector are 
then removed by simply dividing through a heavily Gaussian- 
smoothed version of the mean beam images. This results in 
the desired flats for one observing block, and each COMQ co- 
added frame can be divided by these flats to create calibrated 
on and off beam standard star images, with the pixel-to-pixel 
sensitivity variation removed. This flat fielding procedure was 
carried out separately for each nod position, following which 
the two beams are shifted and averaged (with flux inversion of 
one). 

The target was observed with a large chop angle so that the 
off beam images are devoid of any bright sources. Hence, only 
the flats created from the off beam positions are required to cal- 
ibrate the on beam target images. So, for the target, we simply 
divided each COMQ image by the off beam flat created as de- 
scribed above. This was done identically for the nod position as 
well. Each flat fielded nod position image was then subtracted 
from its immediately preceding flat fielded target image, result- 
ing in 10 calibrated and background- subtracted images of the 
target. One of the point-like sources in each image was used as 

^ http://canadia.ir.isas.jaxa.jp/comics/open/guide/filter 
^ http://www.naoj .org/Observing/DataReduction/ 
^ http://iraf.noao.edu/ 



No.] 



Subaru mid-infrared imaging of the core of M82 



3 



a reference for determining small shifts necessary to create the 
final co-added image. 

Two systematic sources of noise affect COMICS observa- 
tions. Firstly, a low level random noise in the form of hori- 
zontal stripes is introduced into every exposure by the readout 
amplifiers. Although it is largely periodic in nature across the 
field of view as a result of the 16 identical amplifier channels 
employed, this is difficult to remove because of 1) its random 
nature from exposure to exposure; 2) the fact that our target 
covers most of the field of view, leaving free very little of the 
CCD for determination of the noise level. No optimal solution 
was found for removal of this read noise component without 
introducing additional statistical uncertainties in the process, 
so no attempt is made to remove it, except for display purpose. 

Secondly, we find non-negligible background residuals in 
our N11.7 image. This is a result of the fact that 1) the chop 
throw (and nod offset) is relatively large which means that the 
background and target optical light paths differ significantly; 
and 2) the typical time delay between a chopped beam obser- 
vation and its corresponding nod observation is of order several 
minutes, during which time the sky level can vary significantly. 
These residuals are more prominent in the N11.7 image be- 
cause of the wide filter width could make it more sensitive to 
sky variations. In order to remove these, the final Nell image 
was rescaled in order to produce the best match to the target in 
the Nl 1.7 filter, and then subtracted. This leaves a background 
image dominated by the Nl 1 .7 background residuals, with only 
some additional low level structures due to the varying SEDs of 
targets across the field of view between the two wavelengths. 
This background image can then be subtracted from the shifted 
and co-added N11.7 image to create the desired residual-free 
product. 

Final calibrated images are shown in Fig. 1. The attached ab- 
solute astrometry was determined based on multi- wavelength 
comparison of detected sources, as described in the follow- 
ing sections. The coordinates attached to all images presented 
herein are for the J2000 equinox. 

4. Point spread function 

The point- spread-function (PSF) size in our observations 
may be accurately measured from the photometric standard 
star observations, which were carried out both before, and af- 
ter, the target observations. The stellar flux contours are dis- 
played in Fig. 2. Also overplotted is the expected diffrac- 
tion limit resolution for both filters. The 50 per cent stel- 
lar flux contour matches the expected ideal full-width at 
half-maximum (FWHM), meaning that our observations were 
closely diffraction-limited at 0.''36 and in the N11.7 and 
Nell filters, respectively. 

5. Source detection and photometry 

In order to isolate the most significant discrete sources in the 
target field of view, we carried out automatic detection with 
the SExtractor package (Berlin & Arnouts, 1996). The domi- 
nant source of uncertainty for source detection and photometry 
is the diffuse and bright galactic background emission. An it- 



erative process of interpolating the image on a smooth mesh 
and removing detected sources is employed for proper mod- 
eling of this background. But the complex nature of the ob- 
served emission makes the local background quite sensitive to 
the choice of the mesh size. We thus decided to run SExtractor 
in two passes using different mesh sizes. In the first instance, 
a background map is created from a mesh with squares of size 
8x8 pixels. A limit of 3a (with a being the root-mean- square 
variance in the background map at the mesh location of inter- 
est) was chosen for defining significant sources, and a rather 
small minimum detection area of 3 contiguous pixels with a 
significant signal was adopted for detecting the more compact 
sources. In a second pass, we chose a smaller background 
mesh size of 4x4 pixels, and ran the source detection pro- 
cedure again. SExtractor fits ellipse profiles to all sources, 
and only those sources whose centroids remained consistent 
to within some threshold between the two passes were retained 
in the end. A threshold of size equal to the semi-minor axis 
for the Nell filter (and twice this value for the N11.7 images 
which are affected by worse residuals, as described above) was 
found from experimentation to work well. Note that a very 
small background mesh can result in subtraction of some flux 
from the source wings themselves. Test runs on the photo- 
metric star images showed that using a mesh as small as 4 x 4 
pixels results in fluxes lower than the total expected values by 
factors of about 1.2 and 1.3 (N11.7 and Nell, respectively), so 
all SExtractor-determined fluxes have been increased by these 
amounts. Because of the clumpy nature of the galaxy emission, 
additional systematic uncertainties of ^10% are estimated for 
this background correction, and these have also been included. 
A better photometric solution will require accurate mapping of 
the galactic background emission (in longer exposures or using 
space telescopes), which the present data do not allow. 

Calibration of counts to flux was carried out using the pho- 
tometric standard observations described in § 2. Final source 
fluxes are listed in Table 1 and the sources are identified in 
Fig. 3. 

Finally, in order to determine limiting fluxes for non de- 
tections, we simulated two dimensional Gaussian profiles 
matched to the PSF, scaled these according to flux and added 
them in at specific positions in the images where limits were 
required. We then passed these images through our detection 
pipeline, and determined the flux limit corresponding to the 
faintest source that could be detected. We find nominal point 
source flux limits of ^43 mJy and 18 mJy in the Nell and 
N11.7 filters, respectively, at the radio kinematic center posi- 
tion of Weliachew et al. (1984). But these limits can vary sig- 
nificantly with position and source morphology assumed, so 
we quote relevant numbers for some interesting sources indi- 
vidually in the following sections. 

6. Multi-wavelength registration of images 

Registration of images at multiple wavelengths (to better 
than half-arcsec accuracy) is a non-trivial task, given the small 
field-of-view of COMICS, the extended nature of the bulk of 
the emission and the strong dust extinction towards the core. 
Extensive high resolution radio studies over the years have ac- 
curately mapped many supernova remnants (SNRs) and Hll re- 
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Fig. 1. COMICS images in the Nell filter (top) and N11.7 filter (bottom). The field of view displayed is 286x213 pixels wide which, at a pixel scale of 
Of' 133, results in a field size of 38" x28". The linear grayscale is in units of mJy pix~^ . 



No.] 



Subaru mid-infrared imaging of the core of M82 




Fig. 2. PSF measurement from standard stars observed immediately prior to (PI 2UMa), and about one hour after (HD 108381), the target observation. 
The filter central wavelengths (A) are labeled. Contour levels are drawn at 25, 50, 75 and 100 % (last one is a cross) of peak counts. The thick dashed 
gray curve is a circle with diameter l.llXidM (where is the primary mirror diameter of 8.2 m) and is seen to match well with the 50 % contour 
corresponding to the full- width at half maximum. The dotted grid denotes the array pixel scaling of 0'' 133. This shows that the imaging is closely 
diffraction-limited . 
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Fig. 3. Detected discrete sources in the Nell (Top) and N117 (Bottom) filters. Numeric identifications as in Table 1. The source position angle and extents 
are those determined from isophotal photometry with elliptical profiles by SExtractor. 
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gions, which arguably provide the best reference system for 
cross-registration. We used the VLA/MERLIN 5/15 GHz study 
of McDonald et al. (2002, hereafter M02) to search for source 
matches with our mid-IR Nell image. A possible solution was 
identified in which four radio Hll regions coincided with mid- 
IR Nell detections to within Of'5. These four are listed in 
Table 1 and are highlighted in Fig. 4. The rms error of the as- 
trometric solution fit was 0.27 pixels, or Of'036, about 10 times 
smaller than a single resolution element. Adopting this solu- 
tion required a shift of about 0.5 arcsec from the raw Subaru 
astrometry baseline. The astrometry for the N11.7 image was 
tied to that determined for the Nell filter, and all four Hll re- 
gions are also detected in the Nl 1.7 image. 

In order to check whether these mid-IR detections are con- 
sistent with being dust emission from Hll regions, we plot their 
radio to mid-IR spectral energy distributions (SEDs) in Fig. 4, 
and compare to the SEDs of the massive embedded star clusters 
identified by GalHano et al. (2008) in NGC 1365. The overall 
match is excellent given the long 'lever arm' between the radio 
and infrared, and given that some natural variation amongst the 
strengths of the [Nell]Al2.81/im line dominating the Nell fil- 
ter is expected. 

Such a registration procedure is not foolproof, so we have 
carried out cross-checks of our astrometry against published 
(comparatively) low spatial resolution data mentioned previ- 
ously: the [Nell] line map of AL95 and 12.4 jim imaging from 
TG92. A near-perfect astrometric match is found with the for- 
mer at the center of our field of view, but with an increasing ra- 
dial offset towards the edges. The 12.4 fim images of the latter 
show a small systematic offset to larger right ascension with 
respect to both our astrometry and to the line map of AL95. 
All of the above relative offsets lie below ^l'/4 and are un- 
likely to be significant, because the absolute position accura- 
cies or spatial resolution available to these authors were ~1.''5. 
The above comparisons are self-consistent with our astromet- 
ric solution, and we estimate the absolute positional accuracy 
of our images to be better than I". Further detailed checks must 
await high resolution data and more cross-identifications over 
a larger field of view. 

7. Results 

7.7. Comparison with near-IR images 

We start with a qualitative comparison of the overall dis- 
tribution of stars with that of dust. For the stellar distribu- 
tion, we used a 1.6 jim NICMOS image obtained from the 
Hubble Space Telescope archive, and its astrometry was cal- 
ibrated by using the association found by Kong et al. (2007, 
see their Fig. 8) for their point source designated as J095551.0 
in Chandra X-ray imaging. Fig. 5 presents our Nell image and 
the NICMOS 1.6 jim image. These have been overlaid on to 
a single image in Fig. 6. There is a distinct anti-correlation 
between the visibility of stars at near-IR wavelengths, and the 
appearance of mid-IR dust. In particular, the group of well- 
known super-star clusters MGG-7, 9 and 1 1 (cf. Fig. 8 of Kong 
et al. 2007) is neatly nestled within the mid-IR gaps. Other 
clusters, e.g. MGG-3, 6, 8, q and other massive agglomera- 
tion of stars (McCrady et al., 2003), all appear where the mid- 
IR emission is weakest. Heavy and patchy dust extinction is 



known to affect the core of the galaxy - Ay >25 mags (Willner 
et al., 1977; Rieke et al., 1980; O'Connell et al., 1995) - which 
may easily be higher within individual clumps. Using the stan- 
dard interstellar extinction law (Rieke & Lebofsky, 1985), the 
corresponding H band extinction is expected to be Ah>4 mags. 
Thus the effect of dust is sufficient to cause the anti-correlation 
between mid-IR and near-IR structures. 

7.2. Comparison with X-ray images 

We also carried out a comparison against two long X-ray 
images obtained from the Chandra archive: sequences 600735 
and 600736 with observation dates of 2009 Jun 24 and 2009 Jul 
1, respectively. In each case the center of M82 lies on the S3 
chip of Advanced CCD Imaging Spectrometer (ACIS). Results 
presented here use CIAO v4.2 and the CALDB v4.3.0 calibration 
database. The data were re-calibrated using VFAINT cleaning, 
with random pixelization removed and bad pixels masked, fol- 
lowing the software 'threads' from the Chandra X-ray Center 
(CXC)4 to make new level 2 events files. Only events with 
the recommended grades of 0,2,3,4,6 were used. The observa- 
tions were free from large background flares, and after removal 
of time intervals when the background deviated more than 3cr 
from average, the exposure times for sequences 600735 and 
600736 were 118.413 and 118.054 ks, respectively. The de- 
fault astrometry attached to the images was found to agree to 
within <0f'3 with the data presented by Kong et al. (2007), and 
no further refinement was done. In Figs. 5 and 6, an exposure- 
corrected 1.2-5 keV merged image of the two observations is 
included. Note that the brightest source in the field ULX X-1 
is likely to be piled up (see also Feng & Kaaret, 2010) and 
a 'halo' -like counts depression visible in Fig. 5 seems to be 
caused by this. There is also a faint readout streak extending 
from the source position in both directions on the ACIS data. 
But these effects are not relevant for our work, because the chip 
regions used for spectral analysis (later in § 8) do not overlap 
with this; nor do we use this source for astrometric calibration 
or any other analysis. 

There are hints of association of the mid-IR wind stream- 
ers with structures in the diffuse X-ray emission at faint levels. 
There is also a conspicuous dust lane of extinction which bi- 
sects the stellar distribution (e.g. O'Connell et al., 1995), and 
this appears to be associated with cold gas obscuration as well. 
This is manifest as a sharp decrement in the diffuse galactic X- 
ray emission (easily visible in the bottom panel of Fig. 5) co- 
incident with the near-IR deficit. The mid-IR emission avoids 
several spots along this band (in particular, near the dynamical 
center) which may be sites of extremely thick and cold absorb- 
ing material. Detailed analysis of these individual structures is 
left for future work. 

With regard to the X-ray point sources, the most important 
coincidence with an mid-IR source is that of source #18, which 
will be discussed in detail in § 8. The well-known bright ultra- 
luminous X-ray sources (ULXs) known in the core of M82 do 
not show bright mid-IR counterparts. There is, however, one 
potential match. Kong et al. (2007) note that their X-ray source 
J095551.0 (which we use for the NICMOS astrometric caH- 
bration in the previous section) coincides with the radio source 

^ http://cxc.harvard.edu/ciao 
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Fig. 4. Nell (left) and N117 (right) filter zoom-ins of radio Hii regions with assigned IR counterparts listed in Table 1. The four cross-identified sources 
are labelled in red, and their radio positions are marked with I" diameter circles. The corresponding mid-IR counterparts are labelled with brackets in 
green (see Table 1 and Fig. 3. A fifth possible association of the radio source 42.21+59.2 is with mid-IR source # 11; see text in § 6. (Bottom) Radio 
to mid-IR SEDs of the four cross-identified sources are shown. The dashed red and orange lines shows the SEDs for massive clusters M5 and M6 from 
GaUiano et al. (2008). 



42.21+59.2 [B1950] fromM02. Our source # 11 (151.01+45.7 
[J2000]) lies at a separation of Of'5 from the radio position. 
This offset, combined with the fact that the mid-IR counterpart 
is elongated (roughly along the north-south direction; Fig. 3), 
makes cross-identification uncertain. But there are two note- 
worthy points: 1) M02 also found an extension for the radio 
counterpart along a similar direction (see also Fenech et al. 
2008); 2) our mid-IR counterpart bridges the near-IR+X-ray 
position and the radio one. So we tentatively assign this source 
as the fifth Hll region identified in the mid-IR (it is marked in 
Fig. 4). As discussed by Kong et al. (2007), the X-ray counter- 
part may be a ULX hosted within the star cluster detected by 
M02 in the radio (and also now by us in the mid-IR). 

7.3. Discrete sources 

The general source population in the core of the galaxy con- 
sists of discrete sources whose emission is superposed over a 
diffuse radiation field. Of all detected sources, only one (source 
# 8 [150.44+48.2] in Table 1) appears to be truly point-like 
based on the structural parameters returned by SExtractor and a 
FWHM of close to that expected for the diffraction limit. Note 
that the morphology of source #18 [152.70+45.9], which we 
will discuss further in § 8, is uncertain due to its faintness and 
the fact that it is affected by low level read pattern noise. 

Most of the discrete sources appear extended at our high an- 
gular resolution of <0f'4 (corresponding to a linear scale of 
<6-6.5 pc). This is in contrast to the optical and near-IR ap- 



pearance of super star clusters detected in HSTNICMOS imag- 
ing, where most were found to have half light physical extents 
of --3.5 pc or less (O'Connell et al., 1995; McCrady et al., 
2003). Our detected dust features may then be a result of large 
scale outflows from the starburst regions. 

The measured source centroids and fluxes are listed in 
Table 1, and the sources are identified in Fig. 3. There is good 
overall agreement in the emission structures between the two 
images, though the Nell filter shows a greater number of sig- 
nificant individual detections as compared to Nl 1.7, especially 
towards the outer parts of the field of view. Many sources 
have Nell filter fluxes higher than the corresponding N11.7 
filter fluxes by factors of several at least, meaning that the 
[Neii]Al2.81/im dominates in these cases. This is true for all 
the Hll regions cross-matched with the radio catalog of M02, 
which is consistent with the fact that this emission line is a 
strong star-formation indicator, and also why more sources are 
detected in the Nell filter. 

In Table 1 , we comment on some source properties and as- 
sociations with ionized gas sources identified by AL95. In par- 
ticular, we have found additional structure to their 'El, Wl 
and W2' peak emission sites on the eastern and western limbs, 
and on the western ridge of [Nell] line emission. A couple of 
sources are also found on the faint 'bridge' of emission con- 
necting the eastern and western limbs of the ring. Finally, two 
of the tabulated sources (#19 and #21) have Nell fluxes of 
less than 1.5 times the continuum flux in the Nl 1.7 filter, much 
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# 


Mid-IR designation 


^Nell 


^NllV 


Cross-ID 




TO Ann 




mJy 


mJy 


and comments 


1 


148.05+40.5' 


139.30+54.0 


110±12 


19±3 


39.29+54.2 (M02) 


2 


14o.o9+4z.3 


Tin n A , c c n 


95±12 


- 


- 


3 


T /I r\ C . AC o 1" 

149.65+45.2' 


140.91+58.8 


419±44 


80±8 


40.95+58.8 (M02) 


4 


149.68+43.0 


140.93+56.6 


lOlill 


- 


- 


5 


149.92+42.4"^ 


141.18+56.0 


338±36 


85±9 


41.17+56.2 (M02); W2 faint (AL95) 


6 


150.24+43.1 


141.51+56.7 


771±81 


233±24 


W2 bright (AL95) 


7 


150.42+46.2 


141.69+59.8 


212±23 


- 


- 


8 


T C f\ A A AO ^ 

150.44+48.2 


141.70+61.8 


38±4 


27±3 


weak NelliNl 1.7 


9 


150.75+47.3 


142.01+60.9 


72±8 


- 


- 


10 


150.80+44.9"^ 


142.06+58.6 


150±16 


64±7 


42.08+58.4 (M02) 


11 


151.01+45.7 


142.28+59.4 


233±25 


129±13 


Wl faint (AL95); tentative association with 42.21+59.2 (M02) 


12 


151.44+47.3 


142.71+61.0 


321±34 


82±9 


on 'bridge' 


13 


151.44+43.7 


142.72+57.3 


707±75 


153±16 


Wl bright (AL95) 


14 


151.53+49.1 


142.80+62.7 


69±8 


- 


on 'bridge' 


15 


152.11+49.4 


143.38+63.1 


559±59 


145±15 


another component of Wl bright (AL95) 


16 


152.24+47.5 


143.52+61.2 


55±8 


- 


- 


17 


152.50+49.2 


143.78+62.9 


465 ±49 


75±8 


western limb of 'bubble' (El of AL95) 


18 


152.70+45.9 


143.99+59.6 


42±6 


- 


44.01+59.6 (M02); AGN candidate 


19 


TCO Qf\,AC\ Q 

ljz.oU+4y.3 


I44.U0+OJ.U 


29±4 


36±4 


weak NelLNl 1.7 


20 


153.01+47.7 


144.29+61.4 


493±52 


123±13 


southern limb of 'bubble' (El of AL95)) 


21 


153.44+49.6 


144.73+63.3 




31±3 


weak NelLNl 1.7 


22 


154.28+54.2 


145.57+67.9 


142±16 







Table 1. Parameters of detected discrete sources. Source designation is relative to (09h55m, +69d40m) and (09h51m, +69d54m) in J2000 and B1950, 
respectively. The prefix T in the designation refers to the fact that these are Infrared coordinates. The final column states cross-identification with the 
radio catalog of M02, or with diffuse structures in AL95. In some cases, other comments are given. ''^These sources were used for astrometric calibration. 



less than in other cases. This may suggest atypical SEDs with 
continua rising to shorter wavelengths. Upon examining Fig. 3, 
though, it seems that extended emission is present in the Nell 
filter at the positions of both these sources, but SExtractor is 
simply unable to isolate it from the highly clumpy diffuse back- 
ground. In fact, simulating the addition of a weak point source 
at the position of source 21 in the Nell filter suggests a limit of 
~90 mJy which, in turn, implies a NeII/Nll.7 flux ratio con- 
sistent with the distribution of flux ratios for detected sources 
in Table 1 . 

7.4. Diffuse features and galactic emission 

Underlying the discrete sources is a diffuse emission field 
from the galaxy which completely dominates the radiated flux. 
Integrating all observed photon counts over the COMICS field 
of view yields total fluxes of 55 and 108 Jy in the N11.7 and 
Nell filters, respectively. Systematic uncertainties from detec- 
tor cosmetics and the standard star absolute calibration domi- 
nate here (§ 5), for which we allow for 10% variations along 
the full field. The 11.7 fim flux is within a factor of 2 of the 
continuum flux reported by Beirao et al. (2008) for their 'to- 
tal' aperture. This comparison is likely to be full consistent, 
given the difference in aperture sizes and positions used, per- 
haps some contribution of the 11.2 /im polycyclic aromatic hy- 
drocarbon feature to our broad-band photometry, and absolute 
cross-calibration uncertainties. Removing the combined flux 
of the discrete sources identified in Table 1 gives 53 and 106 Jy 
for the two filters. This diffuse emission possesses substantial 
sub-structure, of which we highlight a few specific aspects. 

The most striking features of our high resolution imaging 



are several wind 'streamers', which stand out better in Fig. 7 
which has been smoothed by convolution with a square top hat 
kernel 5x5 pixels (0%1 x 0%1) wide, in order to enhance faint 
diffuse features. Most of these are visible on the northern side, 
with an orientation roughly perpendicular to the major axis of 
the galaxy. The longest of these in our image is a chimney-like 
straight structure, with a linear extent of over 1", or about 120 
pc. Chopping in the mid-IR cuts out a lot of extended emis- 
sion. So, we are undoubtedly seeing only the inner regions of 
much larger streamers and chimneys which have been detected 
on wider scales (e.g. Engelbracht et al., 2006; Kaneda et al., 
2010). 

Towards the north-east of the field of view, two structures 
extend outwards and curve in towards one another, resembling 
the edges of a limb-brightened bubble. Assuming a circular 
shape, the projected diameter of this feature is ^8'', corre- 
sponding to a physical size of ^140 pc. Limb-brightening is 
much less pronounced near the top, which may suggest wind 
rupture, if real. Ruptured bubbles are a natural consequence 
of over-pressurization of outflowing hot gas in the superwind 
model which explains the burst of star formation in the core of 
M82 (e.g. Chevalier & Clegg, 1985; Heckman et al, 1990). If 
there is entrained clumpy dust within this outflow, the dust dis- 
tribution will follow that of the gas. On the other hand, there is 
no obvious evidence from multi- wavelength data that this bub- 
ble is filled in with hot plasma. Its reality as a bubble or as 
two separate streamers need confirmation from further obser- 
vations. 

These various features do not point back to a single site of 
energy ejection. The starburst in M82 is known to occur in a 
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Subaru mid-infrared imaging of the core of M82 




Near-IR 
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Fig. 5. Subaru, HST and Chandra images of the core of M82, aligned 
to match coordinates. The cross is the radio kinematic center deter- 
mined by WeHachew et al. (1984). North is up and East to the left. 
Each image has a field of view of 35''x23''. The annular depression 
of counts visible around the brightest Chandra source to the West is an 
artifact (see 7.2). 



'ring' , rather than a compact nuclear concentration (e.g. Nalcai 
et al., 1987, AL95), and the mid-IR emitting hot dust indeed 
traces out portions of this ring. Multiple burst events of dust 
expulsion along the entire inner ~500 pc region of this ring are 
required to explain the data. 

7.5. Dust temperature and mass 

The mass of the hot, mid-IR emitting matter may be esti- 
mated by assuming optically-thin thermal radiation from uni- 
form spherical dust grains. The largest uncertainty in this es- 
timate is the unknown and spatially- variant dust temperature. 
Telesco & Harper (1980, hereafter TH80) were able to fit a 
modified black body emission model with a temperature of 
~45 K for the dust emitting at far-IR wavelengths of 41-141 
jim over the central '-«^900 pc region. They also found the ob- 
served fluxes at shorter wavelengths from 2-21 jim to overes- 
timate the thermal model predictions by many orders of mag- 
nitude, meaning that hotter dust components characterize the 
mid and near-IR emission. 

In Fig. 8, we overplot the 5-41 jim fluxes and 45 K ther- 
mal model from TH80, along with our integrated N11.7 filter 
flux. There is some non-negligible uncertainty resulting from 
varying beam sizes and centroids amongst the data plotted in 
this figure so we do not carry out a detailed fit. But a quali- 
tative comparison is reasonable if the flux is mostly centrally- 
concentrated. If a two-temperature-only model is assumed, we 
find that a modified black body with flux density varying as 
V^ ^By T (where B^ t is the Planck function at frequency v and 
temperature T) is a fair representation of the 10.5-21 /im inte- 
grated source fluxes for r~160 K. The emissivity power-law 
index of 1.5 has been assumed to be identical to that found for 
the cooler dust by TH80. This model is also plotted in Fig. 8 
(normalized to our 1 1.7 /im data), and shows that the 5 jim ex- 
cess requires other even hotter dust components, if this emis- 
sion is thermal in origin. 

Using this model, the dust mass may be estimated by em- 
ploying the infrared emissivity relation of Hildebrand (1983), 



3 Qv^a TlByj 



(1) 



where Md,D,a and p are the dust mass, galaxy distance, grain 
radius and specific density, respectively. Fy and Qy^a are the 
observed flux density and the grain emissivity. Assuming emis- 
sion from the grain surface into all hemispherical solid angles 
gives the factor of n. Taking canonical values of the size and 
density for a graphite grain composition gives 



for By T expressed in W m~^ Hz~^ sr~^ A mean 11.7 /im 
emissivity of Qy^a^O.009 for graphite grains with diameter of 
^2=0.1 /im is used (Draine, 1985). The observed total flux and 
model temperature then imply Mj~950 Mq for the hot dust. 
In comparison with the mass of cooler 45 K dust determined 
by TH80, this is about 500 times lower. 
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Red: Subaru COMICS (12.81 fjm) 
Green: HSTNICMOS (1.6 fjm) 



Magenta regions : VLA SNRs ( + )&HII(0) 




Fig. 6. RGB overlay of the three panels from Fig. 5, with radio regions plotted as dashed circles (Hii regions) and SNRs (plus signs). The white x sign 
is the kinematic center (Wehachew et al., 1984). There are two heavy plus signs: about 4'' to the East of the center, and 1'' South is the AGN candidate 
in magenta (see § 8), and about 2" to the West and I" South is SN 2008iz in white (Brunthaler et al., 2009). The white box marks the position of the 
unusual radio transient reported by Muxlow et al. (2010). 




54s 53s 52s 51s 9h55m50s 49s 48s 



RA (J2000) 

Fig. 7. Nell filter image smoothed with a convolution filter in order to 
highlight diffuse features. Some of the main features are labelled: 'a' 
and 'b' resemble the projected limbs of a bubble, and 'c' is a long chim- 
ney-like structure. The contours are plotted in a 'square-root' scaling 
in order to highlight faint features, with 10 levels ranging from 1.3 to 
24 mJy pixel" ^ . 
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Fig. 8. Archival wide beam fluxes of the integrated emission from the 
central regions of M82 from Telesco & Harper (1980) (beyond 40 /im) 
and Rieke & Low (1972) (5, 10.5 and 21 /im) overplotted with our 
integrated COMICS 1 1.7 /im flux (red box). The green and blue curves 
represent modified black bodies with temperatures of 45 K (cf. Fig. 1 
of TH80) and 160 K, respectively. The black curve is the sum of these 
two. 
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8. Candidate AGN 

Our target # 18 (152.70+45.9 [J2000]) has a uniquely inter- 
esting story. In our astrometry-corrected images, it is found to 
be coincident (within ^0.''1) with the radio source 44.01+59.6 
[B1950] from McDonald et al. (2001, and references therein). 
Based upon several Hnes of evidence, it has been suggested 
that this may be a weak active galactic nucleus (AGN) in M82. 
The evidence includes the atypical radio SED, detection of OH 
maser satellite lines and also a possible elongated radio jet 
(Wills et al. 1997; Seaquist et al. 1997; Wills et al. 1999). Other 
detailed studies have shown that its SED is not atypical for a 
SNR, and it has an expanding shell with an expansion velocity 
of 2700±400 km s"^ (Allen & Kronberg, 1998; Fenech et al., 
2008). Several hard X-ray AGN searches have been carried out 
over the years (Tsuru et al., 1997; Ptak & Griffiths, 1999), and 
a potential counterpart has been detected with ROSAT (Stevens 
et al., 1999, cf. their source X-3) and even Einstein (Watson 
et al., 1984, who also suggested this as a possible nucleus, 
based on its proximity to the 2.2 /im core). But there is consid- 
erable diffuse soft X-ray emission in this region, and the high 
resolution of Chandra is required for resolving the complex un- 
ambiguously (e.g. Matsumoto et al., 2001). No unambiguous 
and strong X-ray iradio association to our knowledge has been 
reported to date. The source lies close to, but not exactly at 
the dynamical center of the galaxy. It lies 4f'l (70 pc) and 2f'l 
(36 pc) from the radio and optical kinematic centers quoted by 
Weliachew et al. (1984, see also O'Connell & Mangano 1978; 
errors are ^lf'5-2'' in each coordinate). In summary, the na- 
ture of this object still remains uncertain. 

In Fig. 9, we present zoom-in images of the source posi- 
tion in the mid-IR (our Nell filter), the near-IR (NICMOS 
1.6 /im), and in X-rays (Chandra ACIS 0.5-2 keV and 3-7 
keV). There is no point-like source in the near-IR, but it is de- 
tected significantly in our Nell image. A source also clearly 
appears in X-rays, especially at hard energies. We designate 
this source as CXOU J095552.7+694046 (with J2000 coordi- 
nates of 09:55:52.7+69:40:45.8), or J095552.7 for short. Its 
position relative to other X-ray sources within the core of M82 
discussed by Kong et al. (2007) is shown in the bottom panel 
of Fig. 5. We now present analysis of the new long-exposure 
Chandra archival data analyzed herein. In the next section, we 
will also examine the broad-band SED and comment on the 
source nature. 

8. 1 . Chandra X- ray data 

First, we extracted a radially-averaged surface brightness 
profile of the source, and compared this to a simulated PSF im- 
age. The PSF was constructed using the Chandra ray tracing 
package CHART^ assuming observational parameters including 
nominal and source coordinates, roll angle, and ACIS-S chip 
identical to those for the analyzed data (the sequence 600735 
was used for this comparison). The ray trace was then pro- 
jected onto the detector plane and a PSF events file generated 
using MARX. The profile of the source for photons extracted 
over the hard 3-7 keV energy range, compared to the PSF ex- 
pected for 4.5 keV photons, is shown in Fig. 10. A constant 



background level corresponding to the outermost radial bin has 
been removed. The source appears to be consistent with a point 
source, at least within the core. There may be evidence of 
an extended 'wing' around I" radius, but systematic effects of 
the (position-dependent) underlying galactic background can- 
not be ruled out. 

Next, the X-ray spectrum of the source was extracted. A cir- 
cular aperture of radius 1 '.'25 was used to accumulate source 
counts. The diffuse galactic flux around the source is clumpy 
and affected by cold gas absorption to its immediate north. 
An irregularly- shaped polygonal region surrounding the source 
was selected for background subtraction, avoiding nearby point 
sources as well as strong absorption to the north-west. These 
regions are shown in Fig. 9. For spectral extraction the CI AO 
task ps extract was followed by mkacisrmf and mkarf 
to apply the latest calibrations. Models were fitted jointly to 
the two separate observations. Fitting was performed in XSPEC 
(Arnaud, 1996) and includes absorption along the line of sight 
in our Galaxy fixed at a column density of A^h= 4x 10^^ cm~^ 
(based on the data by Dickey & Lockman, 1990)^. 

A minimum grouping of 30 counts per spectral bin was 
applied for fitting and uncertainty determination with the 
statistic. A fit to the resultant 0.7-6.8 keV spectrum (ignoring 
'bad' bins and energies below 0.5 keV) with a single power 
law component results in a hard power law with photon index 
r=l.ll±0.08 (all uncertainties are quoted at 90% confidence), 
but with an unacceptable goodness of fit ;^^=175 for 100 de- 
grees of freedom (dof), with residuals suggesting the presence 
of emission features below 4 keV including those from highly 
ionized Si, Ne and Ar. Several models, including a combina- 
tion of absorbed power laws, or a power law combined with a 
thermal (APEC) plasma model (Smith et al., 2001) were tried. 
Although some of these yield statistically-acceptable fits, the 
power law slope above 3 keV turns out to be very soft, with T 
restricted to ^2-3 in all cases. Furthermore, residuals suggest- 
ing additional emission bands still remain unaccounted for. 

The best model is a combination of two absorbed ther- 
mal plasmas, with temperatures of 0.6 keV and 2.6 keV, ab- 
sorbed by significantly different columns of A^h=8x10^^ and 
4.3 X 10^^ cm~^, respectively. With abundances fixed at Solar, 
assuming the Lodders (2003) abundance table for both temper- 
ature components, a ;^^/dof= 10 1.5/96 is found. The spectrum 
with this best fit is shown in Fig. 10 and the parameters are 
listed in Table 2. Adopting only a single absorber produces an 
unacceptable fit. Regarding abundances, assuming a yield ap- 
propriate for SN II instead (Nomoto et al., 2006, averaged over 
a Salpeter initial mass function from 10 to 50 M©, with a pro- 
genitor metallicity of Z=0.02) produces a fit similar to the best 
one. Letting the abundances of the two components vary also 
does not improve significantly. On the other hand, those for 
SN la (Iwamoto et al., 1999, their W7 model) yield an unac- 
ceptable fit. Residuals suggest a high abundance for specific 
elements, but these are not strictly required by the fit with the 
present data quality. 

The observed flux is Fo.5-10 kev=8.9(±0.4)x 10-^"^ erg s'^ 
cm~^, corrected for Galactic absorption only. Isolating the 
low and the high temperature thermal plasma components. 
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and correcting also for the respective fitted intrinsic obscuring 
columns, yields Fo.5-10 keV=l-2(±0.1)x 10~^^ erg s~^ cm~^ 
and fb.5-10 kev=2.4(±0.2)x 10~^^ erg s~^ cm~^, respectively. 
Though these two values are of the same order of mag- 
nitude, the high temperature component completely domi- 
nates above ~2 keV. Errors refer to those for 90% con- 
fidence intervals determined using the CFLUX convolution 
model in XSPEC. Adding these fluxes gives a total deab- 
sorbed flux of /^o.5-io kev=3.6(±0.2)x 10~^^ erg s~^ cm~^, or 
Z^o.5-iOkeV=5.3(±0.3)xl038 ergs-1. 

Lastly, we point out that the diffuse galactic emission un- 
derlying the source has a soft spectrum, resulting in some ad- 
ditional systematic uncertainty in the source morphology and 
soft band spectral parameters. The source appears to have a 
relatively mild peak excess in Fig. 9. We carried out an addi- 
tional check of the soft band spectrum by changing the back- 
ground extraction regions to an annulus encompassing only 
background counts from the immediate source vicinity. This 
did not remove the need for the soft thermal component, though 
the confidence intervals of the fitted parameters did change, as 
expected. The hard band, however, is not affected significantly 
by the galactic background. 

8.2. Fe K emission line 

The high temperature plasma should also emit a strong ion- 
ized Fe line around 6.7 keV, near the limit of the binned en- 
ergy range available. In order to analyze this line, in Fig. 10 
we show the ungrouped spectrum around 6.7 keV, overplotted 
with the best fit model determined above (binning is applied 
for display purposes only). The default model clearly produces 
a line feature also seen in the data. 

We estimate its significance in a couple of ways. 
Accumulating net source counts over 6.5-6.9 keV after sub- 
tracting off the model continuum measured at the ends of this 
energy range, and comparing these counts with the total (back- 
ground inclusive) counts in the same range, yields a simple 
Poisson signahnoise ratio of 2.4 and 3.1 for the line feature in 
the two Chandra data sets, respectively. Another estimate can 
be obtained by parametrizing the continuum as a power-law lo- 
cally, and then fitting the line as a Gaussian superposed on this. 
The energy range of 5-8 keV is used in order to estimate the 
continuum, and we assume a Gaussian width g fixed at 10 eV. 
The C- statistic is the appropriate one to use for ungrouped data 
(Arnaud, 1996). We find a line energy of 6.68(±0.03) keV and 
a 90% equivalent width range of 0.6-4 keV; in other words, the 
line is significantly detected (despite large uncertainties on its 
strength) and its center energy is consistent with highly ionized 
Fe at 6.7 keV, and not the neutral line energy of 6.4 keV. 

8.3. Comparison with prior archival data 

Detection of any flux variability can be important for un- 
derstanding the nature of the source. The fluxes inferred from 
the two Chandra observations above are consistent with each 
other, and one must examine longer timescales. Searching in 
the Chandra archive for other data, one of the earliest ACIS-S 
observations covering the position of this source is sequence 
600270 with observation date 2002 Jun 18. The exposure time 
after flare removal is only 18.03 ks, but there is a long time 
baseline of about 7 years to the 2009 observations analyzed in 
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Spectral fit for the X-ray counterpart to source #18 



Model components 



Value 



PHABS 




4x10^^ cm-^ 


WABSi 




7.9l};?xl02i cm-2 


APECi 


kTi 
Abundance 
normi 


O.S^tl^l keV 
1 

7.42+1^^x10-^^ 


WABS2 




4.28+}-^^ X 10^2 cm-2 

— 1.05 


APEC2 


Abundance 
norm2 


2.561° j^keV 
1.99t^^xl0-4t 



Table 2. X-ray joint spectral fit to the two Chandra archival data sets 
analyzed. The model is PHABS (wabSi*apeCi + wabS2*apeC2), 
with abundance fixed to that of Lodders (2003). Errors correspond to 
90% parameter uncertainties on a single interesting parameter, when 
not fixed. This model returns a statistic of ;^^/dof=101.5/96. ''^The units 
of norm are 10^^ cm~^. 



previous sections. 

We reduced and analyzed this archival data set in an identical 
manner to the 2009 observations. The extracted net spectrum 
is shown in Fig. 11. Fitting with the same absorbed double 
thermal model as in § 8.1 yields parameters consistent with 
those listed in Table 2 (though with much larger uncertain- 
ties), and a flux Fo.5_io=1.061q2X 10~^^ erg s~^ cm~^ cor- 
rected for Galactic absorption only. The 90% confidence in- 
terval overlaps with the flux determined from the 2009 obser- 
vations (§ 8.1). 

Thus, there is no indication for strong flux variability on 
a timescale of several years either. A similar inference has 
also been reached by Chiang & Kong (2011) in a recent 
analysis of other observations (their source designation is 
CXOU J095552.8+694045). Our comparisons limit any flux 
variations to less than 20%. 

9. Discussion 

We have presented the highest resolution imaging of the nu- 
clear regions of M82 to date in the mid-IR. The Subaru diffrac- 
tion limit at wavelengths of 11.7 and 12.81 jim is 0.''36 and 
Of' 39, respectively. Compared to the previous best resolution 
observations (the 12.4 {im imaging with PSF of If'l presented 
by TG92), our images have PSFs improved by factors of 3.1 
and 2.8 at 1 1.7 and 12.81 /im, respectively, and we cover a core 
region larger by a factor of at least two. The referee also made 
us aware of a conference proceeding by Ashby et al. (1994), 
where an image at 11.7 /im with an angular resolution of 0% 
is presented. The absolute astrometry of their image agrees to 
within ~1'' with our work, as well as with the works of AL95 
and TG92. To our knowledge, no further details have been 
published. 

Multi-wavelength image comparisons show an anti- 
correlation between the observed stellar distribution (probed in 
the near-IR) with the distribution of warm dust that we probe. 
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Fig. 9. Magnified images (9"wide x 6" high) on the position of the AGN candidate (source # 18) in the mid-IR (top left; COMICS Nell filter) and near-IR 
(top right; HST NICMOS F160W filter). The X-ray images at the bottom are for energy ranges of 0.5-2 keV (left) and 3-7 keV (right), respectively. 
Because the PSF of ACTS is larger than in the IR, both X-ray images are zoomed out by a factor of two; the dashed box in the bottom-right panel shows 
the field size of the IR panels. The scale bar notches denote 1" offsets from the source position in both RA and Dec. The regions used for spectral 
extraction of the target and background are shown in the 0.5-2 keV image as the magenta circular and polygonal apertures, respectively. A secondary 
annular background region used for a consistency check is also shown in yellow. 



This means that obscuration in these dusty regions is heavy 
enough to scatter even near-IR radiation. 

Our observations provide the best view of the base of the 
dusty superwind that is known to exist in this galaxy, and reveal 
several elongated features on projected physical scales of up to 
£=120 pc at least. If the mid-IR emitting dust is mixed in with 
and entrained in outflowing gas, then the travel time (t) is 

i \ /200kms-i^ 



^ = 5.9 



X 10^ yr. 



(3) 



120 pc 

where a velocity of v=200 km s~^ identical to that found by 
Nakai et al. (1987) for the molecular gas is used as reference. 
This time period suggests recent energy input from young star- 
bursts. The kinetic energy required to expel this dust is only a 
fraction of that channeled into the total gas mass, and so expul- 
sion by SNe occurring at various locations around the ring of 
star-formation can easily account for this (Nakai et al., 1987). 

Using the integrated fluxes over the COMICS field-of-view 
alongside archival data at longer and shorter wavelengths, and 
assuming only a two temperature phase for the dust, we are 
able to describe the broad-band mid-IR diffuse emission with 
a r=160 K modified black body with a hot dust mass of 
~1000 M0, in addition to the 45 K cool dust component 
responsible for the far-IR emission. Assuming a standard 
gas: dust ratio of 100, the mass of the gas associated with the 
mid-IR emitting dust is then ~10^ Mq. This is lower than the 
ionized gas mass of 2 x 10^ M© (TH80; Willner et al. 1977) by 
a factor of about 20, meaning that the ionized gas may be the 
predominant environment for the hot dust that we are observ- 
ing. The cooler dust, on the other hand, is likely to be tracing 
the distribution of molecular gas (TH80). 



More than 20 discrete sources are detected, and most are 
found to have extended profiles. Matching against radio cat- 
alogs suggests that we have resolved at least four (and tenta- 
tively, five) Hll regions in the mid-IR for the first time. These 
Hll regions have monochromatic continuum dust luminosities 
ranging from vL^'^ ^"'=2x 10% (radio ID 39.29+54.2) to 
8x lO^L© (41.17+56.2), consistent with being powered by em- 
bedded super star clusters similar to those seen in some other 
nearby galaxies (e.g. Galliano et al., 2008). 

As seen in Fig. 6, no source is detected at the position of 
the radio transient identified by Brunthaler et al. (2009) as 
SN2008iz, the brightest radio SNR in M82 over the past 20 
years. The start of bright radio flaring activity is limited to the 
time interval of 2007 Oct 29 to 2008 Mar 24. Our observa- 
tions were carried out one month after the final set of follow- 
up radio observations on 2009 Apr 04 reported by Brunthaler 
et al., when they found the source to have an integrated 22.2 
GHz flux of 9.2±0.2 mJy. Our flux limits are FNeii<38 mJy 
and Fnii.7<18 mJy, assuming a point source of angular size 
equal to the diffraction limits in each filter. In general, there 
is little overlap between confirmed SNRs and mid-IR detec- 
tions. Additionally, no mid-IR counterpart is detected at the 
location of the unusual radio transient found by Muxlow et al. 
(2010). The source appeared on radio images obtained within 
the period of 1-5 May 2009 - contemporaneous with our mid- 
IR observations - and was not present one week before. It is 
located ^1 from the position of our source #18 (see also 
Kong & Chiang 2009 for an identification of the X-ray coun- 
terpart). We estimate point source upper limits of 34 mJy and 
15 mJy in the Nell and Nl 1.7 filters. 
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Fig. 10. s (Top) Azimuthally-averaged Chandra ACIS-S radial surface brightness profile of J095552.7 (points with error bars). The empty gray squares 
show the normalized counts profile before background subtraction. The filled black circles denote the normalized net profile, and these can be compared 
to the expected PSF from a point source (histogram). The source profile is extracted from the 3-7 keV image, and the PSF from ray-tracing of 4.5 keV 
photons. (Bottom Left) X-ray spectrum and residuals to a fit with two absorbed thermal plasma models. Data from both analyzed observations (black 
and red for sequences 600735 and 600736, respectively) are included. The two thermal models are shown as the dotted lines. (Bottom Right) Ungrouped 
spectrum around the Fe K line and best fit model from the left panel (rebinned to obtain a minimum 2a significance for plotting only). 



No space observatory is foreseen to have a better resolv- 
ing power than Subaru at ~10 /im, though the Mid InfraRed 
Instrument (MIRI) onboard JWST should provide excellent 
sensitivity for only a modest loss in angular resolution at the 
same wavelengths (e.g. Wright et al., 2004). MIRI is also 
expected to have a field of view larger by a factor of ^2 
on a side, making source identification much more secure. 
On a longer timescale, the European/JAXA mission Spica 
(Nakagawa, 2004) will be crucial for deep searches of far-IR 
signatures of AGN activity at this source location. On the 
ground, the Gran Telecopio Canarias, with a primary mirror 
diameter of 10.4 m, can improve the resolution slightly to Of'3 
under good observing conditions. Even better resolution must 
await larger ground-based observatories such as the Extremely 
Large Telescope. This will not be diffraction-limited under 
natural seeing conditions and will require additional adaptive 
optics capabilities in the mid-IR to improve upon the results 
presented herein. 



9. 1 . Nature of AGN candidate source 

The puzzling radio source 44.01+59.6 has a significant Nell 
mid-IR detection (our source #18 [152.70+45.9]), and is also 
coincident with a source visible in high-resolution X-ray im- 
ages (J095552.7), to well within our estimated absolute po- 
sitional uncertainties. From Chandra data, it is immediately 
apparent that the X-ray spectrum is not consistent with that 
of AGN, which usually display broad-band power-laws with 
photon-indices characteristic of radiatively-efficient ac- 
cretion (e.g. Mateos et al., 2005). Low luminosity AGN with 
radiatively-inefficient flows or jets, too, usually display hard 
X-ray continua (e.g. Yu et al., 2010). Similarly, one can argue 
against direct association with other kinds of accreting sources 
such as X-ray binaries. 

The detection of a strong He-like Fe line instead of a neutral 
line centered on 6.4 keV implies the absence of cold reflecting 
matter, also atypical for an accreting source. Highly ionized 
lines have actually been attributed to heavily obscured AGN in 
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Fig. 11. ACIS-S spectrum of J095552.7, extracted from an 18 ks expo- 
sure (sequence 600270) observed on 2002 Jun 18, seven years prior to 
the data shown in Fig. 10, and fitted with the absorbed double thermal 
model described in § 8.1. A minimum bin grouping of 20 counts has 
been applied. 

some cases (e.g. Iwasawa et al., 2005; Nandra & Iwasawa, 
2007, and references therein), though on larger scales. In such 
a scenario, the AGN itself is not readily apparent; rather, the 
visible spectrum is dominated by surrounding gas which is ir- 
radiated by the AGN along directions out of our line of sight. If 
this is true for the case of source J095552.7, the underlying ac- 
creting object (be it an AGN or an X-ray binary) could power 
the radio jet found by Wills et al. (1999) and also provide a 
ready source where accreting gas is in rotation (as inferred from 
the maser lines found by Seaquist et al. 1997). But this would 
also result in X-ray photoionization spectral features or even an 
ionized reflection continuum. An extra layer of optically-thick 
obscuration is then also required (in addition to the two that we 
detect) with an extreme covering factor very close to unity so as 
to completely hide the reflection continuum and any neutral Fe 
K line. This is not consistent with the detection of a jet which 
would be expected to decrease the covering factor by clearing 
away some surrounding matter. 

Finally, the 2-10 keV luminosity of 5x 10^^ erg s~^ means 
that the source is unlikely to have a bolometric X-ray power 
which would place it in the ULX regime (e.g. Makishima et al. 
2000) because of the steep spectra of the fitted thermal models. 

Thus, our analysis allows us to reach some firm conclusions 
as to what the source is not. The true source nature still remains 
uncertain, though, and in the following, we discuss plausible 
alternatives. 

In addition to our [Nell] image, the source appears most 
prominently as a compact object in hard X-rays. One possi- 
bility may then be that both the mid-IR and the hard X-ray 
components are associated with the SNR visible in the radio. 
To test this, we may compare the source power to that of Gas 
A, a powerful and young Galactic SNR. From Suzaku obser- 
vations summed over the entire remnant, the 4-10 keV X-ray 
luminosity of Gas A is determined to be L^^f^^^y^4xl0^^ 
erg s~^ (Maeda et al., 2009). This is about 150 times smaller 
than the deabsorbed X-ray luminosity L4-10 keV^6x 10^^ erg 
s~^ that we find for our source by using the high temperature 
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APEC component alone (§ 8.1). This requires that the electron 
density (Ue) in the X-ray emitting interstellar medium (ISM) 
around source # 1 8 be much larger than in the case of Gas A. rie 
may be determined by using the normalization values returned 
by the APEC2 component and the emission volume which, in 
the scenario under consideration here, should correspond to the 
volume of the radio- emitting plasma. Fenech et al. (2008, see 
their Table 3) measure a diameter of 0.86 pc for the SNR, us- 
ing which yields an electron density of ^e~1.75(±0.3)x 10^ 
cm~^. A low plasma filling factor would only push this value 
up. This is, indeed, much larger than typical ISM densities of 
1-10 cm~^ relevant for Galactic SNRs. We also note that the 
present limit on any X-ray variability in the source (^20% over 
seven years) derived from archival data comparisons in § 8.3 is 
consistent with the observed smaller flux changes in Gas A over 
similar timescales (e.g. Patnaude et al., 2010). The detected 
coUimated radio jet may then suggest some kind of axisymme- 
try in the progenitor (e.g. a binary merger event). 

On the other hand, the mid-IR power for source # 18 
(AL;l (12.81 /im)=l-5>< 10^^ ^^g s~^) is in excess of the integrated 
12 iim IRAS power of Gas A (^L^^fif ^m)=5x 10^^ erg s-^; 
Saken et al. 1992) by a much larger factor of ^3000. But 
given the unknown dust temperature, the fraction of freshly 
formed vs. ambient heated dust, and the fact that our Nell fil- 
ter flux is likely to be dominated not by continuum but by the 
[Nell]Al2.81/im line (which is known to be a strong cooling 
line for SNRs; e.g. Rho et al. 2008), further comparison of the 
mid-IR fluxes is difficult. 

The other possibility to consider for the X-ray thermal 
plasma components is that of a hot ISM phase in a compact star 
cluster. In fact, collisionally ionized plasmas with a wide range 
of temperatures have been previously inferred to exist within 
starburst galaxies (e.g. Iwasawa et al., 2005; Ranalli et al., 
2008; Strickland & Heckman, 2009), though on larger scales. 
Fig. 12 shows the compiled radio-to-X-ray fluxes for source # 
18. A mid-IR flux dominating the SED is consistent with an 
origin in a starburst, though this is not a unique solution. The 
detection of two distinct layers of absorption (§ 8.1) may also 
suggest a scenario which combines the above possibilities as 
follows. The hot X-ray thermal plasma and the radio counter- 
part could both be associated with a SNR, which is embedded 
within a host star cluster. The high column density affecting 
the APEC2 component (Table 2) may easily be explained by 
strong local absorption within molecular clouds, for instance. 
The point-like profile of the hard X-ray data is also consistent 
with this. The cluster, on the other hand, could appear promi- 
nently in soft X-rays and also in the mid-IR. The comparatively 
low aborption affecting this component is then due to gas along 
the plane of the galaxy on larger scales. 

Further insight into the nature of this source will be possi- 
ble if it can be isolated at longer wavelengths characterizing 
the peak of typical star-formation SEDs (this may be within 
reach of JWST), or through detection of spectral features in 
the sub-mm with ALMA. Subsequent long Chandra exposures 
would be useful to place tighter constraints on X-ray variabil- 
ity. Meanwhile, the question of whether M82 hosts an AGN or 
not remains to be answered. 
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Fig. 12. Broad-band SED of source #18. The radio data are the 1.3 
cm to 74 cm fluxes (and one limit) from Allen & Kronberg (1998). In 
the mid-IR, we plot our Nell filter flux, and a 3a Nl 1.7 detection limit 
measured assuming Poisson statistics within an aperture equal in size 
to the diffraction limit. The X-ray regime shows the unfolded spectrum 
relevant for the model fitted to the observed data in Fig. 10. 
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